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Results from the NA49 experiment at the CERN SPS are presented on event-by-event transverse 
momentum and multiplicity fluctuations of charged particles, produced at forward rapidities in 
central Pb-|-Pb interactions at beam momenta 20A, 30A, 40A, 80A, and 158A GeV/c, as well as 
in systems of different size {p-\-p, C-|-C, Si-|-Si, and Pb-|-Pb) at 158A GeV/c. This publication 
extends the previous NA49 measurements of the strongly intensive measure ^pj, by a study of 
the recently proposed strongly intensive measures of fluctuations A[Pt,N] and E[Pt,V]. In the 
explored kinematic region transverse momentum and multiplicity fluctuations show no significant 
energy dependence in the SPS energy range. However, a remarkable system size dependence is 
observed for both A[Pt, V] and S[Pt, V], with the largest values measured in peripheral Pb-|-Pb 
interactions. The results are compared with NA61/SHINE measurements in p -|-p collisions, as well 
as with predictions of the UrQMD and EPOS models. 
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PACS numbers: 25.75.-q, 25.75.Gz 


I. INTRODUCTION AND MOTIVATION 


Ultra-relativistic heavy ion collisions are studied mainly to understand the properties of strongly interacting matter 
under extreme conditions of high energy densities when the creation of the quark-gluon plasma (QGP) is expected. 
The results obtained in a broad collision energy range by experiments at the Super Proton Synchrotron (SPS) at 
CERN, the Relativistic Heavy Ion Collider (RHIC) at BNL, and at the Large Hadron Collider (LHC) at CERN 
indeed suggest that in collisions of heavy nuclei such a state with sub-hadronic degrees of freedom appears when the 
system is sufficiently hot and dense. 

The phase diagram of strongly interacting matter is most often presented in terms of temperature (T) and bary- 
ochemical potential {(J-b), which reflects net-baryon density. It is commonly believed that for large values oi fj,B the 
phase transition is of the first order and turns into a rapid but continuous transition (cross-over) for low ^jlb values. 
A critical point of second order (CP) separates these two regions. The phase diagram can be scanned by varying the 
energy and the size of the colliding nuclei and the CP is believed to cause a maximum of fluctuations in the measured 
final state particles. More specifically, the CP is expected to lead not only to non-Poissonian distributions of event 
quantities like multiplicities or average transverse momentum BE]: but also to intermittent behavior of low-mass 
7r“''7r“ pair and proton production with power-law exponents calculable in QCD Hi]. 

The NA49 experiment at the CERN SPS [S] pioneered the exploration of the phase diagram by an energy scan for 
central Pb-|-Pb collisions in the range 20A to I58A CeV {^/snn = 6.3-17.3 CeV), as well as a system size scan at the 
top SPS energy of 158A CeV. Evidence was found IS El that quark/gluon deconfinement sets in at a beam energy of 
about 30A CeV. Thus the SPS energy range is a region where the CP could be located. At present the search for the 
critical point is vigorously pursued by the NA61/SHINE collaboration at the SPS [HI and by the beam energy scan 
program BES at RHIC [5]. 

The NA49 experiment already measured multiplicity fluctuations in terms of the scaled variance u of the distribution 
of event multiplicity N ESin] and event-by-event fluctuations of the transverse momentum of the particles employing 
the strongly intensive measure [I1E3]. The present paper reports a continuation of this NA49 study by analyzing 
two new strongly intensive measures of event-by-event transverse momentum and multiplicity fluctuations, A[Pt,N] 
and T,[Pt,N] [CTI15j . These measures are dimensionless and have scales given by two reference values, namely they 
are equal to zero in case of no fluctuations and one in case of independent particle production. Unlike they allow 
to classify the strength of fluctuations on a common scale. 

This paper is organized as follows. In Sec. |h] the new strongly intensive measures of fluctuations A[P 7 ’,A^] and 
T,[Pt,N] are introduced and briefly discussed. Data sets, acceptance used for this analysis, detector effects, and 
systematic uncertainty estimates are discussed in Sec. |HI| The NA49 results on the energy and system size dependences 
of transverse momentum and multiplicity fluctuations quantified by the new measures are presented and discussed in 


Sec. IV A summary closes the paper. 


II. STRONGLY INTENSIVE MEASURES OF TRANSVERSE MOMENTUM AND MULTIPLICITY 

FLUCTUATIONS 

In thermodynamics extensive quantities are those which are proportional to the system volume. Examples of 
extensive quantities in this case are the mean multiplicity or the variance of the multiplicity distribution. In contrast, 
intensive quantities are defined such that they do not depend on the volume of the system. It was shown m that 
the ratio of two extensive quantities is an intensive quantity, and therefore, the ratio of mean multiplicities, as well 
as the commonly used scaled variance of the distribution of the multiplicity N, w[Y] = ((iV^) — (7V)^)/(V), are 
intensive measures. Finally, one can define a class of strongly intensive quantities which depend neither on the volume 
of the system nor on the volume fluctuations within the event ensemble. Such quantities can be truly attractive when 
studying heavy ion collisions, where the volume of the produced matter cannot be fixed and may change significantly 
from one event to another. Examples of strongly intensive quantities are mean multiplicity ratios, the $ measure of 
fluctuations [16], and the recently introduced A and E measures of fluctuations [141 [15]. In fact, it was shown [14] 
that there are at least two families of strongly intensive measures: A and S. The previously introduced measure 4> is 
a member of the E-family. 

In nucleus-nucleus collisions the volume is expected to vary from event to event and these changes are impossible 
to eliminate fully. Thus, the strongly intensive quantities allow, at least partly, to overcome the problem of volume 
fluctuations. Generally, the A and E measures can be calculated for any two extensive quantities A and B. In this 
paper B is taken to be the aecepted particle multiplicity, N {B = TV), and A the sum of their transverse momenta Pt 
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(A = Pt = 'Y^^=iPTii the summation runs over the transverse momenta pt; of all accepted particles in a given event). 
Following Refs. [3[T5] the quantities A^] and T,[Pt,N] are defined as: 


- {PtMN]] 


( 1 ) 


and 


S[PT,iV] = + {PtMN] - 2{{PtN) - (Pr)(iV))], 


( 2 ) 


where: 


and 


w[Pr] 


{pM - {Prf 
(Pt) 


w[A] = 


(N^) - (Nf 

(N) 


( 3 ) 

( 4 ) 


are the scaled variances of the two fluctuating extensive quantities Pt and N, respectively. The brackets (...) represent 
averaging over events. The quantity uj{pt) is the scaled variance of the inclusive pt distribution (all accepted particles 
and events are used): 


uj{pt) 


Pt - Pt'^ 

Pt 


( 5 ) 


Equations 0 and 0 can be used only when assuming that oj{pt) is not equal to zero. There is an important 
difference between the A[Pt,N] and "SIPtM] measures. Only the first two moments: (Pt), (N), and {Pt"^), {N'^) 
are required to calculate A[Pt, A], whereas E[Pt, A] includes also the correlation term {PtN) — {Pt){N). Therefore 
A[Pt, A] and S[Pt, A] can be sensitive to various physics effects in different ways. In Ref. [ll] all strongly intensive 
quantities containing the correlation term are named the E family, whereas those based only on mean values and 
variances the A family. As already mentioned, the previously studied [HIS] measure belongs to the E family 
and obeys the relation: 


‘J’pT = \/prw(pT)[\/E[PT, A] — 1] (6) 

With the normalization of A and E proposed in Ref. m these quantities are dimensionless and have a common scale 
making possible a quantitative comparison of fluctuations of different, in general dimensional, extensive quantities. 
The basic properties of the A[Pt, A] and E[Pt, A] measures are the following: 

1. Absence of fluctuations. In the absence of event-by-event fluctuations (A = const.^ Pt = const.) the values of 
A[Pt, A] and E[Pt, A] are equal to zero. 

2. Independent Partiele Model (IPM). If the system consists of particles that are emitted independently from each 
other (no inter-particle correlations) A[Pt, A] and E[Pt, A] are equal to one. For this case vanishes. 

3. Model of Independent Sources (MIS). When particles are emitted by a number (As) of identical sources, which are 
independent of each other and P(As) is the distribution of this number, then A[Pt, A](As) and E[Pt, A](As) 
are independent of As (intensive measures) and of its distribution P(As) (strongly intensive measures). The 

measure has the same property. An example of MIS is the Wounded Nucleon Model (WNM) [T7], where 
As = Nw (number of wounded nucleons). Another example is a model where nucleus-nucleus {A-\-A) collisions 
are an incoherent superposition of many independent nucleon-nucleon (A -|- A) interactions. For these cases all 
three fluctuation measures, namely A[Pt, A] and E[Pt, A], are independent of the number of sources (and 
therefore insensitive to the centrality of the collisions) and have the same values for A + A and A-|-A collisions. 
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The measures A[P 7 ',-^] and T,[Pt,N] have similar advantages to w[A]. Like a;[A] they have two reference values, 
namely w[-/V] equals zero when the multiplicity is constant from event to event and one for a Poisson multiplicity 
distribution. Therefore one can judge whether fluctuations are large ( > 1 ) or small ( < 1 ) compared to independent 
particle production. However, a;[iV] is not a strongly intensive quantity, and in the MIS one finds CLi[iV](As sources) = 
a;[iV] (1 source) + (n)a;[fVs], where (n) is the mean multiplicity of particles from a single source and a;[iV 5 ] represents 
fluctuations of Ns- 

A comparison of the properties of A[Px,N], E[Pr, A], and is presented in Tabled 

The quantities A[Pt,N] and T,[Pt,N] were studied in several models. The results of^imulations of the IPM, the 
MIS, source-by-source temperature fluctuations (example of MIS), event-by-event (global) temperature fluctuations, 
and anti-correlation between Pt/N and N were studied in Ref. [I8]. Predictions from the UrQMD model on the 
system size and on the energy dependence of AlP^jN] and T,[Pt,N] are shown in Ref. [TS]. Finally, the effects 
of quantum statistics were discussed in Ref. [T5]. The general conclusion is that A[Pt,N] and E[Pt,A^] measure 
deviations from the superposition model in different ways. Therefore, the interpretation of the experimental results 
may benefit from a simultaneous measurement of both quantities. 



Unit 

No fluctuations 

Independent Particle Model 

Model of Independent Sources 

^’PT 

MeV/c 

4’pt = -yJpTio{pT) 

4>pr = 0 

independent of Ns and P{Ns)', 
^pj.{Ns sources) = 'l?pr(l source) 

A[PT,iV] 

dimensionless 

o 

II 

<1 

A[PT,iV] = 1 

independent of Ns and P{Ns)', 
A[Pt,N]{Ns sources) = A[Pt,A']( 1 source) 

E[PT,iV] 

dimensionless 

E[PT,iV] = 0 

S[PT,iV] = 1 

independent of Ns and P{Ns)', 
T,[Pt,N]{Ns sources) = E[Pt, 771(1 source) 


TABLE I: Properties of 4>py, A[Pr,A^], and E[Pt,A] in the absence of fluctuations {N = const., Pt = const.), in the 
Independent Particle Model and in the Model of Independent Sources. 


III. DATA SELECTION AND ANALYSIS 

The data used for the analysis, event and particle selection criteria, uncertainty estimates and corrections are 
described in the previous publications of NA49 [HIS] on the measure Here we only recall the key points. 

The analysis of the energy dependence of transverse momentum and multiplicity fluctuations uses samples of Pb-|-Pb 
collisions at 20A, 30A, 40A, 80A and 158A GeV/c beam momenta (center of mass energies from 6.3 to 17.3 GeV per 
N + N pair) for which the 7.2% most central reactions were selected. The analysis of the system size dependence is 
based on samples of p + p, semi-central C-l-G, semi-central Si-|-Si, and minimum bias and central Pb-|-Pb collisions 
at 158A GeV/c beam momentum. Minimum bias Pb-t-Pb events were divided into six centrality bins (see Ref. [T^] 
for details) but due to a trigger bias the most peripheral bin (6) is not used in the current analysis. For each bin of 
centrality the mean number of wounded nucleons {N\y) was determined by use of the Glauber model and the VENUS 
event generator [5D] (see Ref. m)- 

Tracks were restricted to the transverse momentum region 0.005 < pr < 1.5 GeV/c. For the study of energy 
dependence the forward rapidity range 1.1 < y* < 2.6 was selected, where y* is the particle rapidity calculated in the 
center-of-mass reference system. For the study of system size dependence at 158A GeV/c the rapidity was calculated 
in the laboratory reference system and restricted to the region 4.0 < < 5.5 |12j (it approximately corresponds to 

1.1 < y* < 2.6). As track-by-track identification was not applied, the rapidities were calculated assuming the pion 
mass for all particles. For the energy scan an additional cut on the rapidity y* calculated with the proton mass was 
applied (y* < yheam ~0-5) [131 [21]. This excludes the projectile rapidity domain where particles may be contaminated 
by e.g. elastically scattered or diffractively produced protons. 

The acceptance of azimuthal angle (f> was chosen differently for the study of energy and system size dependence 
(Fig.g. For the energy scan a common region of azimuthal angle was selected for all five energies (only particles 
within the solid curves in Fig. (left) were retained), whereas a wider range was used at 158A GeV/c for the system 
size study (see Fig. (right)). Together with the track quality criteria and rapidity cuts this results in using only 
about 5 % respectively 20 % of all charged particles produced in the reactions. 

An additive correction for the limited two track resolution of the detector was applied to the values of A[Pt,N] 
and ^[Pt, A^]- The procedure to determine this correction was analogous to the one used to estimate the corrections 
for in Refs. [HHSj. Mixed events were prepared for each of the analyzed data sets and then processed by the 
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FIG. 1: (Color online) Examples of NA49 {(pjPr) acceptances of charged particles with the azimuthal angle of negatively 
charged particles reflected (see Ref. [TS] for details). The solid lines represent the analytical parametrization of acceptance used 
for further analysis. Left: acceptance used for the energy scan of pt and N fluctuations, example for 2.0 < Pn < 2.2. Right: 
acceptance used for the system size dependence of pt and N fluctuations, examples for 1.2 < < 1.4 and 2.4 < yt < 2.6. 

Additional cut on y* (see the text) not included. Figure reproduced from Refs. | 121113 | where parametrizations of the curves 
can be found. 
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FIG. 2: (Color online) Additive corrections dA[PT, N] (left) and dTi[PT,N] (right) for limited two track resolution in the 7.2% 
most central Pb+Pb events at 20A-158A GeV/c. Estimates for positively charged, negatively charged and all charged particles 
are distinguished by different markers (see legend). 


NA49 simulation software. The resulting simulated raw data were reconstructed and the measures A[Pt,N] and 
T,[Pt,N] were calculated using the same selection cuts as used for the real events. The additive two track resolution 
corrections 6A[Pt, fV] and (511 [P-r, N]) were calculated as the difference between the values of A[Pr, TV] (or S[Pr, .^j) 
after detector simulation and reconstruction and before this procedure. The resulting corrections for the data of the 
energy scan are plotted in Fig. and those for the data of the system size study in Fig. 

The statistical uncertainties on A[Pt,N] and were obtained via the sub-sample method [T^llIS]- The 

systematic uncertainties were estimated by varying event and track cut parameters (the procedures were identical to 
those applied for in Refs. nang). 


IV. RESULTS AND DISCUSSION 

The results shown in this section refer to accepted particles, i.e., particles that are accepted by the detector and 
pass all kinematic cuts and track selection criteria as discussed in Sec. [ml The data cover a broad range in pt 
(0.005 < Pt < 1.5 GeV/c). The rapidity was restricted to the interval 1.1 < y* < 2.6 (forward rapidity) where 
contamination from beam produced 5-rays is small. The selected azimuthal angle region is large and represents 
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FIG. 3: (Color online) Additive corrections 5A[Pt, (left) and 5 E[Pt, N] (right) for limited two track resolution versus mean 
multiplicity (N) of accepted particles for p + p, C+C, Si+Si and three different centralities of Pb+Pb events at 158A GeV/c. 
Estimates for positively charged, negatively charged and all charged particles are distinguished by different markers. The lines 
represent the analytical parametrization: SA{'E)[Pt, N]{{N)) = —a^(N) with a being the parameter of a fit to the data points. 


essentially the whole detector acceptance for the study of the system size dependence at 158A GeV/c (see lines in 
Fig. 0 (right)). It is more limited for the analysis of the energy dependence since the same region was chosen at all 
energies (see lines in Fig. (left)). Results are not corrected for limited kinematic acceptance. Such a correction is not 
possible since it depends on the, in general, unknown correlation mechanism. Instead the limited acceptance should 
be taken into account in the model calculations when comparing to experimental results. However, corrections for 
limited two track resolution of the NA49 detector were applied (see Sec. Ill and Refs. [Tain]). A possible bias due 
to particle reconstruction losses and contamination in the accepted kinematic region was estimated to be small and 
is included in the systematic uncertainty of the results. 


A. Energy scan for central Pb+Pb interactions 

Figure [^presents for the 7.2% most central Pb+Pb interactions the energy dependence of the fluctuation measures 
A[Pt, a] and S[Pt, A] calculated separately for all charged, negatively charged, and positively charged particles. The 
sample of negatively charged particles is composed mainly of tt~ mesons, whereas the sample of positively charged 
particles is dominated by 7r+ mesons and protons. Therefore, the measured values of A[Pt,N] or T,[Pt,N] could 
differ between both charges. Moreover, among all charged particles additional sources of correlations could exist that 
are not present in positively or negatively charged particles separately. For all three charge selections the values of 
A[P 7 ’, a] are smaller than one, the expectation for independent particle production. For T,[Pt, A] fluctuations for all 
and positively charged particles are close to the hypothesis of independent particle production (similar to the results 
on m which belongs to the same family of strongly intensive measures), whereas for negatively charged particles 
S[Pt, A] values are higher than one. It was suggested in Refs. [HI [19] that values of A[Pr, A] < 1 and S[Pr, A] > 1 
can be explained as due to effects of Bose-Einstein statistics. Similarly, > 0 was predicted in Refs. dH [53] as a 
consequence of Bose-Einstein correlations. 

The measured values of A[Pt, A] and E[Pt, A] are compared to predictions of the UrQMD PH |2S] and EPOS PH 
[27] models in Fig. (solid and dashed lines respectively). The models do not simulate a phase transition or the 
critical point. However, resonance decays and effects of correlated particle production due to energy-momentum, 
charge and strangeness conservation laws are taken into account. The most central 7.2% interactions were selected for 
comparison of the energy scan results, in accordance with the real NA49 events. The procedure of selecting the 7.2% 
most central events was the following: a sample of minimum bias Pb+Pb events was produced. Then the distribution 
of the impact parameter b was drawn and the value of b was determined below which 7.2% of the events remained. 
The resulting impact parameter range was 0 < 6 < 4.35 fm in UrQMD and 0 < 6 < 4.00 fm in EPOS. Finally, high 
statistics samples of UrQMD and EPOS events were produced in these impact parameter ranges separately for each 
energy. 
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FIG. 4: (Color online) Energy dependence of A[Pt, N] (top) and S[F’t, N] (bottom) for the 7.2% most central Pb+Pb interac¬ 
tions. Statistical uncertainties are denoted by lines, systematic ones by color boxes. Data (points) are compared to predictions 
of the UrQMD 3.4 (solid lines) and EPOS 1.99 (dashed lines) models with acceptance restrictions as for the data. 


The measures A[Pt,N] and E[Pt,N] were calculated from charged particles, consistent with originating from the 
main vertex. This means that mostly pions, protons, kaons and their anti-particles from the primary interaction 
were used because particles coming from the decays of Kg, A, E, S, O, etc. are suppressed by the track selection 
cuts. Therefore, the analyses of UrQMD and EPOS events were also carried out by using primary charged pions, 
protons, and kaons and their anti-particles. The tracking time parameter in the UrQMD model was set to 100 fm/c 
and therefore the list of generated kaons, pions and (anti-)protons did not contain the products of weak decays. In 
the parameter settings of the EPOS model the decays of Kg^j^, A, E, S, O, etc. particles were explicitly forbidden. 
Einally, in the analysis of the UrQMD and EPOS events the same kinematic restrictions were applied as for the NA49 
data. 

Eigure (top) shows that the energy dependence of A[Pt,A^] in the UrQMD model exhibits behavior similar to 
that observed in the measurements. In both cases one finds A[Pt,N] < 1, i.e. values below those for independent 
particle production. As Bose-Einstein correlations are not implemented in the UrQMD model we conclude that in this 
model there must be another source(s) of correlation(s) leading to A[Pt,N] < 1. The EPOS model shows A[P 7 ’, A] 
values which are significantly higher that those obtained from the NA49 data and UrQMD. The comparisons for 
T,[Pt, Af] can be seen in Fig. (bottom). Here the predictions of UrQMD he above the measurements for all charged 
and positively charged particles, whereas they are significantly below the results for negatively charged particles. On 
the other hand EPOS calculations for negatively charged particles are close to the data, but exceed the measurements 
even more than the UrQMD predictions for all charged and positively charged particles. 

The measured energy dependences of A[Pt, A] and AiIPt, A] do not show any anomalies which might be attributed 
to approaching the phase boundary or the critical point. However, it should be noted that due to the limited acceptance 
of NA49 and the additional restrictions used for this analysis the sensitivity for such fluctuations may be small if the 
underlying range of correlations in momentum space is large. 


B. System size dependence at 158A GeV/c 


Figure [^presents the dependence of A[Pr, A] and E[Pr, A] at 158A GeV/c on the size of the colliding nuclei as 
well as on the centrality of Pb-|-Pb interactions. The measured values for all accepted charged particles and also for 
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FIG. 5: (Color online) Dependence of A[Pt, Af] (top) and E[Pt,A'^] (bottom) on the mean number of wounded nucleons (Nw) 
on the size of the colliding nuclei {p - triangles, C - stars, Si - squares, Pb - circles) and the centrality of Pb+Pb interactions 
at 158A GeV/c. Statistical uncertainties are denoted by error bars, systematic uncertainties by colored boxes. Data (points) 
are compared to predictions of the UrQMD 3.4 (solid lines) and EPOS 1.99 (dashed lines) models with acceptance restrictions 
as for the data. 


positively and negatively charged particles separately are plotted versus the mean number of wounded nucleons. The 
values of A[P 7 ’, A?'] for p + p, C+C, Si+Si, and the two most central classes of Pb+Pb collisions are lower than one. 
For the more peripheral Pb+Pb interactions A[Pt,N] increases above one to a maximum for the most peripheral 
Pb+Pb collisions. The values of 'S[Pt,N] for negatively and all charged particles are significantly above unity (the 
prediction of the independent particle production model) and also reach a maximum in the most peripheral Pb+Pb 
interactions. For positively charged particles the values are close to zero or below. The same behavior was observed 
for the measure [12] ^[Pt,N] belong to the same family of strongly intensive measures). Finally, it 

is worth recalling that also for multiplicity fluctuations a maximum was observed in peripheral Pb+Pb collisions by 
NA49 [m 

Figure]^ shows that A[PT,fV] and E[Pt, A^] for all charged particles are usually higher than for either negatively 
or positively charged particles. Moreover, in case of T,[Pt,N], the values for positively charged particles are always 
lower than those for the negatively charged particles. In Ref. [T2| it was shown that also values of for positively 
charged particles were lower than those for negatively charged and for all charged particles. However, the same effect 
was observed in simulations using the HIJING model and the fact that values for positively charged particles were 
always lower than those for negatively charged ones was found to be related to the limited acceptance and treatment 
of protons as pions in the calculation of rapidity. 

To further investigate the nature of the correlations leading to the observed values of the measures A[Pt,N] and 
Ti[Pt,N] a toy model was constructed, in which an anti-correlation between mean transverse momentum per event 
{Pt/N) versus multiplicity {N) was assumed [18]. The parametrization of Pp/N versus N was taken from the NA49 
p + p data [121 (in the current paper the same p + p data are used), resulting in A[Pt, A^j = 0.816(0.005) and 
S[Pt, A^] = 1.008(0.002). This shows that Bose-Einstein correlations are not the only candidate for the explanation 
of A[P 7 ', A^] < 1 and E[Pt, A^] > 1 [HI [12], but this observation, especially for smaller systems, may also be explained 
as due to the known Pp /N versus N anti-correlation [18] . 

The system size dependences of A[Pp,N] and T,[Pp,N] were also compared to predictions of the UrQMD and 
EPOS models (the procedure of selecting the proper impact parameter range was analogous to that used in the case 
of the energy scan). Figure!^ shows that in the UrQMD model the values of A[P'r, A^] and S[Pt, A^] are only weakly 
dependent on {Nw}- For ZTfPTjA^] the values are slightly below independent particle production and for S[Pt,A^] 
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they are close to 1. The pronounced maximum of fluctuations seen in the data at (Nw) ~ 90 in Pb+Pb collisions 
is not reproduced by the UrQMD model. The EPOS model predictions for A[Pt,N] are similar to the predictions 
of UrQMD. Also the trends in I1[Pt,N] for negatively charged and positively charged particles are quite similar in 
UrQMD and EPOS. In contrast, for all charged particles the values of 'E,[Pt,N] are much higher in EPOS than in 
UrQMD and describe the NA49 results surprisingly well. 


C. Search for the critical point 

When searching for possible indications of a critical point it is most appropriate to plot the strength of fluctuations 
using the standard phase diagram coordinates temperature T and baryochemical potential fiB- Moreover, central 
collisions of nuclei provide the cleanest interaction geometry. Eor such reactions fits of the hadron gas model (see 
e.g. Ref. [55]) were performed to determine the temperature Tchem and baryochemical potential /is of the produced 
particle composition. These values are believed to be close to those of the hadronization along the transition line in 
the phase diagram. The value of Tchem was found to decrease somewhat for collisions of larger nuclei, whereas 
decreases rapidly with collision energy. 

Results for A[Pt,N] and T,[Pt,N] for inelastic p + p as well as central Pb+Pb collisions are shown in Fig. [^versus 
Pb- The p + p results from NA61 [29ll30], plotted for comparison, were obtained using the NA49 acceptance cuts. 
One observes little dependence on pB for both Pb+Pb or p + p collisions. In particular, there is no indication of a 
maximum that might be attributed to the critical point. A similar conclusion was reached from the pB dependence 
of US]. The measurements of A[Pt,N] are consistent for the two reactions. The values of 'E[Pt,N] are close to 
unity with the exception of the higher result in Pb+Pb for negatively charged particles. 

The dependence of A[P 7 ’,iV], and T,[Px,N] on Tchem is shown in Fig. [7] at the beam momentum of 158A GeV/c 
for p + p, semi-central C+C, Si+Si and central Pb+Pb reactions. The results ior p + p from NA49 (solid triangles) 
and NA61 (open triangles) are consistent. A maximum is observed for Si+Si interactions similar to the one found 
previously for ^p^ in Ref. [I2]. There it was interpreted as a possible effect of the critical point (31] consistent with 
QCD-based predictions of Ref. mm- Interestingly, for the same system, studies of intermittency in the production 
of low mass 7r“'"7r“ pairs |33| and of protons [34] found indications of power-law behavior with exponents that were 
consistent with QCD predictions for a CP. 

Unfortunately, theoretical predictions, for fluctuations at CP, are not yet published for the new fluctuation measures 
A[P 7 ’, N] and ^[Pt, A^]- However, calculations for Si+Si collisions at 158A GeV/c by using the Critical Monte Carlo 
(CMC) model [31 |3S] are currently under study. 


V. SUMMARY 

This paper reports on the continuing search at the CERN SPS by the NA49 experiment for evidence of the critical 
point of strongly interacting matter expected as a maximum of fluctuations. Results are presented on transverse 
momentum and multiplicity fluctuations of charged particles, produced at forward rapidities (1.1 < y* < 2.6) in 
central Pb+Pb interactions at beam momenta 2QA, 30A, 40A, 80A, and 158A GeV/c, as well as in different systems 
{p + p, C+C, Si+Si, and Pb+Pb) at 158A GeV/c. New strongly intensive measures of fluctuations, A[Pt,N] and 
Ti[Pt,N], were measured. This paper is an extension of previous NA49 studies [T2l|T3] where the strongly intensive 
measure was used to determine transverse momentum fluctuations. The quantities A[Pj-,N] and S[P 7 ',iV] are 
dimensionless and have two reference values, namely they are equal to zero in case of no fluctuations {Pt = const., 
N = const.) and one in case of independent particle production. Therefore, A[Pr, A^] and TiIPtjN] are preferable 
to 4)p,j, for which only one reference value is defined, i.e. ^pj. = 0 MeV/c for the model of independent particle 
production (IPM). 

The NA49 results show no indications of a maximum in the energy dependence of transverse momentum (see also 
Ref. [21]) and previously measured multiplicity [3l] fluctuations in central Pb+Pb collisions throughout the CERN 
SPS energy range (but finer steps in the scan would be desirable). At all energies the values of A[Pt,N] are below 
one, i.e. below the expectation from the IPM. The values of T,[Pt,N] are close to one for all charged and positively 
charged particles and slightly higher in case of negatively charged particles. The effect for negatively charged particles 
can probably be explained as due to Bose-Einstein statistics. For positively charged particles the interpretation is less 
clear because other sources of correlations (for example resonance decays) can contribute to the correlation measures. 

The system size dependence of A[Pt,N] and T,[Pt,N] (and the related measure [H] 131]) at 158A GeV/c 
shows maxima for Si+Si and peripheral Pb+Pb interactions (for A[Pt, A] this maximum is mostly seen for peripheral 
Pb+Pb interactions). In central collisions values of A[Pt, A"] are lower than one, whereas values of T,[Pt,N] for all 
charged and negatively charged particles are higher than one (weak anticorrelations). The maximum of T,[Pt, A^] for 
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FIG. 6: (Color online) Energy (/rs) dependence of A[Pt,A^] (top) and E[Pt,A'^] (bottom) for the 7.2% most central Pb+Pb 
interactions and a comparison to NA61 inelastic p + p interactions. pB values for Pb+Pb collisions are taken from Ref. |28| . 
NA49 data indicate that at the top SPS energy pb does not depend on the system size (C+C, Si+Si, Pb+Pb). Therefore, 
the pb values for p + p are also displayed and assumed to be the same as for Pb+Pb. NA61 data were taken from Refs. |29II30| . 
For NA61 only statistical uncertainties are shown. 


all charged particles in C+C and Si+Si interactions is about 5% higher than the base line defined by the IPM. Also 
previously studied multiplicity fluctuations for the most central A + A collisions were found to show a maximum for 
Si+Si reactions at 158A GeV/c [31) . The excess of transverse momentum and multiplicity fluctuations is two times 
higher for all charged than for negatively charged particles as expected for the CP [T]. 

The NA49 collaboration also searched for evidence of the critical point in an intermittency analysis of low-mass 
7r+7r“ pair [33] and proton [31| production. Indications of power-law behavior consistent with that predicted for a 
CP were found in the same Si+Si interactions at 158A GeV/c. The intriguing results strongly motivate the ongoing 
critical point search by the successor experiment NA61/SHINE |5] which performs a systematic two-dimensional scan 
(SPS energies and system size {p, Be, Ar, Xe, Pb)) of the phase diagram of strongly interacting matter. A maximum 
of several GP signatures, the so-called hill of fluctuations, would signal the existence of the CP. The RHIC Beam 
Energy Scan pursues a complementary program measuring higher order moments and cumulants of net-charge 
and net-proton distributions in Au+Au collisions. So far no clear evidence for the CP was found [331137) . Thus the 
possible existence of the CP remains an interesting and challenging question. 
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